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Role of gluons and the quark sea in the proton spin * 

Petr Zavad43 

Institute of Physics AS CR, Na Slovance 2, CZ-182 21 Prague 8, Czech Republic 

The real, interacting elementary particle always consits of a ’bare’ particle and a cloud of virtual 
particles mediating a self-interaction and/or the bond inside a composite object. In this letter we 
discuss the question of spin content of the virtual cloud in two different cases: electron and quark. 
Further, the quark spin is discussed in the context of proton spin, which is generated by the interplay 
of quarks and virtual gluons. We present a general constraint on the gluon contribution and make 
a comparison with the experimental data. 

PACS numbers: 12.39.-x ll.10.Gh 13.60.-r 13.88.+e 


1. INTRODUCTION 

In our recent paper 1J we studied the proton spin structure in leading order of the covariant approach assuming the 
gluon contribution to the proton spin can be neglected. However the question of the real role of gluons in generating 
the proton spin is still open. Actually some recent data obtained at RHIC and their analyses nn can suggest a 
positive gluon contribution to the proton spin. The main aim of the present letter is to extend discussion from our 
previous study to the case of nonzero gluon contribution. We will show what constraint on the gluon contribution 
follows from the covariant approach. 

Sec. [2]is devoted to a discussion about some general aspects of particle spin and its scale dependence. Two different 
examples are considered, electron and quark. The electron spin structure is also an interesting topic, see the recent 
study [7j and the previous papers 019]. Particularly important questions concern the spin of quarks inside the 
nucleon. In Sec. [3] the discussion about the proton spin, which is generated by the interplay of angular moments of 
quarks and gluons, continues in the context of recent experimental data. 

In present calculations we use, as before, the rest frame of the composite system as a starting reference frame. This 
frame is suitable for the consistent composition of spins and OAMs of the constituents in the representation of spinor 
spherical harmonics. The resulting state serves as an input for construction of the covariant quantities, like the spin 
vector or the spin structure functions [T]. 


2. SPIN OF THE PARTICLE IN ITS SCALE DEPENDENT PICTURE 

In general, description of real interacting particles can be related to their ’bare’ or ’dressed’ form. In our present 
discussion we address the general questions: 

a) How much do the virtual particles surrounding bare particle contribute to the spin of corresponding real, dressed 
particle? 

b) How much do the virtual particles mediating bond of the constituents of a composite particle contribute to its 
spin? 

In quantum mechanic the total angular momentum (AM) of any particle including composite ones is given by the 
sum of the orbital AM (OAM) and spin, J = L + S. The corresponding quantum numbers are discrete sets of integers 
or half-integers and in the relativistic case only total AM conserves, so only J and J z can be the good quantum 
numbers. We will illustrate the problem with two different examples, electron and quark. 


2.1. Spin of electron 


The electron, as a Dirac particle, in its rest frame has AM defined by its spin, s = 1/2. This value is the same 
for the dressed electron (as proved experimentally) and for the bare one (as defined by the QED Lagrangian). The 
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FIG. 1: Scale dependent image of a real particle, see text. 


dressed electron is a bare electron surrounded by the virtual cloud of 7 and e _ e + pairs, as symbolically sketched 
in Fig. [l] for different scales represented by the parameter Q 2 . So the renormalization as a continuous change of 
the scale should not change the AM represented by the discrete numbers, J e (Q 2 ) = s = 1/2. But what about the 
projections J|(Q 2 )+ J](Q 2 ) = ± 1 / 2 ? Can the contribution of virtual cloud J](Q 2 ) differ from zero and how much? 
In this letter we present a semiclassical estimate of the vector J 7 . The electromagnetic field, or its 7 -quanta, are 
according to Maxwell equations created by the electric current. We consider the current generated by the electron 
states represented by the spinor spherical harmonics 


1 3,3*) = ( r ) 


1 

V2e 
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where w represents the polar and azimuthal angles ( 0 , tp) of the space coordinates r with respect to the axis of 
quantization z, l p = j ± 1 / 2 , X p = 2 j — l p ( l p defines the parity), energy e = y/k 2 + m 2 and 
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where l p = j — 1/2 and X p = j + 1/2. The functions Yi i_ are usual spherical harmonics. The radial functions in the 
case of free electron read: 


Rki{r) = \J‘ 27 ^J l+1/2 (kr), (3) 

J r 2 RkiRk'idr = 2ix5{k - k'), 

where k = |k| and J u (z) are Bessel functions of the first kind, otherwise (e.g. for electron in the hydrogen atom) 
the radial functions differ according to an external field. However, it is important that the information about the 
electron AM is completely absorbed in the angular terms and does not depend on the radial functions. The states 
| j,j z ) are eigenstates of the total AM and have been discussed before [T] in momentum representation, while in the 
present note we deal with their coordinate representation corresponding to the rest frame of an composite system. 
The corresponding current reads 


= ( J o,I) = &] lpja (r) 7°7( r ) 


(4) 


and one can check that 


I 0 = hipj >jz (cos 9) , I = hnp jdz (cos 9) r, 


(5) 
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3—3 cos 20 

4 

3 ± 1 

2 ’ 2 

5+3 cos 29 

4 

5 1 5 

2 ' 2 
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64 
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64 
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TABLE I: The examples of the angular distributions Pj,j z - The common factor 1/An is omitted. 
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( 6 ) 


A few examples of pjj z are given in Table [TJ where one can observe the following. The stationary current + depends 
only on j and \j z \ } therefore it does not involve any information on the direction of electron polarization. So, there 
is no reason to expect any correlation between electron polarization and polarization of the electromagnetic field 
generated by this current, or equivalently polarization of the statistical set of emitted and reabsorbed 7 . In other 
words the average polarization of virtual cloud of 7 and consequently also e _ e + pairs should be zero. The AM of the 
electromagnetic field is given by the relation 


J 7 = J r x (E x H)d 3 r, (7) 

where E, H are the corresponding intensities of electric and magnetic field. Due to the symmetry of current ([5]) that 
generates these fields, the corresponding AM satisfies 


J 7 = 0, ( 8 ) 

the proof is given in Appendix. This relation follows only from the angular terms in the wave function |T]) and does 
not depend on the radial ones. In other words, the relation ([ 8 ]) holds not only for a free electron, but also for an 
electron bound in the hydrogen atom. The result represents a mean value, which is not influenced by the fluctuations 
generated by single 7 . So, this calculation suggests the integral AM of the cloud of virtual 7 is zero despite the fact 
that AM of its source, the electron in a state ([l]), is not zero. While the free electron emits and reabsorbs virtual 
photons by itself, the electron bounded in hydrogen atom in addition exchanges (emits and absorbs) virtual photons 
with the proton. Since the AM of the electromagnetic field generated by the proton is zero as well, the total AM 
of hydrogen will be given only by AMs of the electron and proton, without contribution of the electromagnetic field 
generated by both the particles. 

Similar arguments can be relevant also for atoms in general and perhaps for the nucleons bound in a nucleus. This 
would suggest the virtual particles mediating the binding of nucleons also do not contribute to the resulting spin of 
nucleus, which must be always integer or half-integer. 

Our approach (A) has a common basis with the QED calculation (B) suggested in Ref. 0, since both the approaches 
follow from the general QED equations (1),(6) and (7) in the last reference. Despite that, there are some differences 
between them, the most apparent are as follows. 

(A) The approach is semiclassical only and the AMs are directly related to the electron wavefunction ip and the 
classical electromagnetic field A M generated by the electric current ip'y^ip. The preferred reference frame is the frame 
defined by the spinor spherical harmonics ([!]) or the rest frame of the defined composite system (e.g. atom). This 
simplified treatment allows us to obtain the relation ([ 8 ]) but without a decomposition into the spin and OAM parts. 

(B) The study is focused on the fundamental problem of the AM decomposition in quantum field theory and explicit 
calculation is performed for the QED. The preferred reference frame is the infinite momentum frame. The light-front 
formalism is adopted to achieve a compatibility with the standard formulation of the parton model. The result on 
the total boson (photon) AM seems be also rather small, Sb + Lb ~ <D(e 2 ). 

So both the approaches, adopting rather different formalism and correspondingly also some different assumptions, 
are not contradictory. From a phenomenological point of view, they are complementary. 
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FIG. 2: Scale dependent image of a valence quark, see text. 


2.2. Spin of quark 

The situation with quarks inside a nucleon is more complicated. The quark at different scales is sketched in Fig. [l] 
(in which the bare electron surrounded by virtual cloud of 7 and e - e + pairs is now replaced by the bare quark with 
a cloud of virtual g and qq pairs). The terminology is as follows: 

i) The bare quark can be identified with the current quark, which can be described by the distribution functions 
q a (x ) defined in the quark-parton model. They are related to the sets of quarks and antiquarks in the figure for 
Q 2 —> 00. 

ii) The constituent quark can be identified with the dressed quark at a low Q 2 scale. 

iii) The valence quark can be identified with the set of quarks, from which the cloud of virtual qq pairs and gluons 
is separated off. In the figure the valence quarks are represented by the central spots. Strictly speaking, depending 
on the scale, valence and sea quarks may not be clearly distinguishable. In a short time interval At, a quark from 
the virtual qq pair is indistinguishable from the source, valence quark, see Fig. [2j However, the usual definition in 
terms of the quark-parton model distributions 


<ivai( x ,Q 2 ) = <f( x ,Q 2 ) - q a {x,Q 2 ) (9) 

is unambiguous. For quarks the parameter Q 2 represents the renormalization scale, but also the DIS parameter 
(— Q 2 = photon four-momentum square) or equivalently a scale of the space-time domain inside which the photon 
absorption takes place [T]. 

Now, we can put the question b) for the quarks bound inside the proton: How much the field of virtual gluons 
generated by the valence quarks and the sea of virtual qq pairs created by the gluons contribute to the proton spin? 
This question is being studied in the experiments, which measure contribution of the gluons and sea quarks to the 
proton spin. Available data from the experiments COMPASS [TO] and HERMES m suggest rather small gluon 
contribution, in fact the data are consistent with zero within statistical errors. The very recent results of the RHIC 
experiments EHB] suggest a positive gluon contribution which, however still cannot fully compensate for a small 
quark contribution to the proton spin. 

These results can be interpreted in the framework of covariant approach presented in Ref. [I], In the paper we 
studied the relativistic interplay between the quark spins and OAMs, which collectively contribute to the proton spin. 
The simplest scenario assuming 

1) the quarks are in the state j = 1/2, see Eq. 113[1], 

2 ) mass of quarks can be neglected, (m/e) -7 0 , 

3) there is no gluon contribution, i.e. proton spin J = 1/2 is generated only by the AM of quarks, see Eq. 109[1: 
gave a prediction for the contribution of the quark spins in DIS region, 

A£=l, ( 10 ) 

while the “missing” part of the proton spin is fully compensated by the quark OAM. This prediction fits the data 
mm surprisingly well. 

However, in a more general case, if only condition 1) is assumed, then AM of each quark consists of the spin and 
OAM part 


1 + 2/t . _ 2 — 2 p, . (l z ) 2 — 2/1 

9 J z 1 \ Z / 9 3 Z1 / \ 1 1 o — 

3 3 (s z ) 1 + 2/i 


( 11 ) 


where j z = ±1/2, see Eqs. (17), (22) in Ref. [1]. It is very important result, which easily follows from the algebra 
of spinor spherical harmonics representing the solutions of Dirac equation. The ratio g, = (m/e) here plays a crucial 
role, since it controls a ’’contraction” of the spin component, which is compensated by the OAM. It is a quantum 
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FIG. 3: Dependence of AS on J g and fl. The dotted, full and dashed lines correspond to AS = 0.25, 0.33 and 0.4 respectively. 

mechanical effect of relativistic kinematics. Mathematically, a small fl means the lower component of Dirac spinor is 
important. This is also the case in mm, where with some distinction of formalism and paradigm the authors come 
to similar results on the spin of quark bound in the proton. The effect is strongly correlated with the transversal 
motion of quarks inside the nucleon m- 

The relations (11) imply that the total AM of a composite system of quarks with j i = j? = j '3 = ... = 1/2 reads 


J q = (S z > + (L z ), 

where the ratio of the total spin (S z ) and OAM (L z ) is the same as for the one-quark states above: 

(L z ) 2 - 2 p 

(§ z ) ~ 1 + 2 p 

Further, if the proton spin consists of the quark and gluon contributions, one can write 




J ' = ^ 


J 9 = \(l-x). 


With the use of Eqs.(12)-(14) one gets 


<S*> + (L z ) = (S z ) 1 + 


2 — 2 p 

1 + 2 fl 


2 *’ 


= 1-2J 9 , 


which after replacing (S z ) = AS/2 gives 


AS = - (1 - 2 J 9 ) (1 + 2/1). 
o 


( 12 ) 


(13) 


(14) 


(15) 


(16) 


In a higher approximation, if we admit an admixture of the quark states with j > 3/2, the relation is modified as 


AS < - (1 - 2J 9 ) (1 + 2fi ). 

o 


(17) 


3. DISCUSSION AND CONCLUSION 


The relation (16 1 means the quark spin content depends on two parameters, the gluon contribution J 9 and the 


quark effective mass ratio fl, which affects proportion of the quark OAM. It follows from the algebra of spinor spherical 
harmonics and from general rules of AM composition in the system of quarks and gluons with the total spin J = 1/2. 
The dependence is demonstrated in Fig. [3j One can observe: 

a) AS <1/3 corresponds to J 9 > 0 for any 1 > jl > 0. A special case AS = 1/3 and fl —0 implies J 9 -A 0. 

b) AS > 1/3, then the sign of J 9 depends on fi. Apparently, J 9 < 0 would imply J q > 1/2. 

In this way the COMPASS and HERMES data n~2TTDT| giving AS « 1/3, can be compatible also with a positive 
gluon contribution J 9 suggested by the recent data on RHIC. Note that positive J 9 correlates with a positive quark 
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effective mass ratio ft. Can we somehow estimate this parameter? In the proton rest frame a quark has momentum 
(k) = -y/3/2 (kx) and there are independent ways to estimate it, for example: 

1) If we take proton diameter d„ = 0.84fm, then the uncertainty relation gives for the corresponding momentum 
roughly k » 230MeV. 

2) The analysis [T5H2T1 of the data on the azimuthal asymmetry in semi-inclusive DIS suggest (kx) ~ 400 — 600MeV. 

3) The statistical approaches f22h24l suggests ( k ) « 40 — lOOMeV and a similar value was obtained in Ref. [25] for 
valence quarks, see also discussion in Ref. EM- 

Therefore, if one suppose the quark effective mass of the order MeV, then the parameter fl < 0.1, which gives a 
similar upper limit on the gluon contribution, J 9 < 0.1. Finally, the spin contribution of the sea quarks is known to 
be small or compatible with zero |14j . This can confirm the expectation that the sea quark contribution correlates 
with the gluon contribution. 

In the present approach the quark effective mass ratio fl and the gluon AM contribution J 9 are free, phenomeno¬ 
logical parameters constrained by the relation (161. At the same time corresponding scale-dependent parameters, for 
example the quark effective masses, are (at least in principle) calculable in QCD [28]. However, due to nonpertur- 
bative aspect of related task, the real calculation can be extremely difficult. That is why our approach based on the 
covariant quark-parton model can be a useful supplement to the exact but more complicated theory of the nucleon 
spin structure based on pure QCD. 

To conclude, interpretation of the available sets of experimental data in framework of the covariant approach 
suggests an important role of the quark OAM for the creation of the proton spin on the scale Q 2 defined by the 
data. The data of experiments COMPASS and HERMES on the quark and gluon contribution to the proton spin 
are fully compatible with our approach. At the same time, a limited positive gluon contribution does not contradict 
the covariant approach. However, the precise data on J 9 are still missing, so the existing experimental data do not 
disprove the hypothesis J 9 ss 0 based on the analogy with AM of virtual photons given by Eq. Ill- 
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Appendix A: Proof of the relation (J8]) 


The current ([5| generates the electric and magnetic field 

E(r) = 

H( r )= / I(r ' )X 


If we define 



h x (r')pj,j z (cos 9') r' 


" 3,0 

|r - r' 


,3/2 


dV, 


S(r) 


hi(r')Pj,j z (cos 9') j3 _, 

, ,3/2 a r 

r — r'| ' 


where X = J, II, then with the use of (§,© we get: 

E(r) = -W J (r) + S(r)r, H(r) = W n (r) x r. 


In terms of spherical coordinates 


(Al) 

(A2) 


(A3) 


(A4) 


ri = r sin 9 cos ip , r 2 = r sin 9 sin ip, = r cos 9 


(A5) 
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we have 


S( r) = 

W A >) = 


hi(r')pjj z (cos 6') 


(r 2 + r' 2 — 2 rr' (sin 0 sin O' cos (<p — ip') + cos 0 cos 0')') 
hx(r')pj,j, (cos O') r' 


(r 2 + r ’ 2 — 2 rr' (sin 8 sin O' cos (ip — p') + cos 0 cos 0')) 
Obviously S(r) does not depend on ip so we have 

S(r) = S(r,0). 

In the second integral, after substitution ip = ip' — ip we replace correspondingly in r' : 

x' = r' sin O' cos ip' —> r' sin O' (cos ip cos ip — sin ip sin tp ), 
y' = r' sin O' sin <p' — ► r' sin O' (cos ip sin ip + sin ip cos ip) 


Yj^r' sin O' dip'dO'dr 1 , 
jj^r' sin O'dip'dO'dr'. 


and instead of (A7) we obtain 


W*<„. [ {«**«»*- ^^mAs*****, 

J ( r 2 _(_ r / 2 _ 2 rr' (sin 0 sin O' cos ip + cos 0 cos 0')) 

W 2 X ( r) = J 

W x (r) = J 


hx (t*)P i i z (cos O') r' sin O' (cos ip sin p “l - sin ip cos p) a . 

-—- YFT r sm " dydO dr , 

(r 2 + r' 2 — 2 rr' (sin 0 sin O' cos ip + cos 0 cos O')') 

h x ( r ')pj,j z (cos O') r' cos O' ,2 . 

- --- j— r sin OdipdOdr. 

(r 2 + r' 2 — 2 rr' (sin 0 sin O' cos ip + cos 0 cos 0')) 


Since in general 


feven(lp) sin Ip dip = 0 , 


(A6) 

(A7) 

(AS) 

(A9) 

(A10) 

(All) 

(A12) 

(A13) 

(A14) 


where f e ven(' l P ) = feven(-ip), then the second term in ( |A11 1 ,( A12) vanishes and the expressions are simplified as 

W?(v) = W x (r,0)n, W x (r) = W x (r,0)r 2 , W x (r) = W x (r, 0), (A15) 

where 


W x (r,0) = 


h-x(r')pj,j z (cos0') r' sin0' cosip 


■ sin 0 


w x (r,o) = j 


(r 2 + r' 2 — 2 rr’ (sin 0 sin 0' cos ip + cos 0 cos 0')) 


Yj^r’ sin O’ dip dO 1 dr 1 , 


hx(r')pj,j z (cos O') r' cos O' 


r' sin O' dip dO'dr'. 


(r 2 + r' 2 — 2rr' (sin 0 sin O' cos ip + cos 0 cos 0')) 3 ^ 2 
After inserting from ( |A4[ ) into ([j]) we integrate the AM density 

j 7 = r x ((—W 7 (r) + S(r)r) x (W 7/ (r) x r)) , 
which with the use of ( |A 8 [ ), ( |A15[ ) gives 

(W I W 11 (r\r 2 r 3 + r 2 r 3 ) + W II W 3 r 2 r 3 — W I W 3 I (r\r 2 + r 2 ) — W 3 W 3 T r 2 r 3 ) , 
j 7 = ^ (— W I W 11 (r±r 2 r 3 + rfr 3 ) — W II W 3 r\r 3 + W I W 3 I (r\r 2 + rf) + W^W^rirs) , 

0 

(-W n r 2 r 3 + M / 3 // r 2 ) , 

+ Sr 2 <| (W H nr 3 - W 3 r r-[) , 

0 


(A16) 

(A17) 

(A18) 

(A19) 


These terms depend on <p only via coordinates ry and r 2 (A5). Since each term involves just one odd power of 
7 "i ~ cost/? or r 2 ~ sint/?, the corresponding integral satisfies (18]). 
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